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Type II quasars are luminous Active Galactic Nuclei whose centers are ob- 
scured by large amounts of gas and dust. In this paper we present 3-band HST 
images of nine type II quasars with redshifts 0.2 < z < 0.4 selected from the 
^ ! Sloan Digital Sky Survey based on their emission line properties. The intrinsic 

luminosities of these AGNs are estimated to be —24 > Mb > —26, but optical 
obscuration allows their host galaxies to be studied unencumbered by bright nu- 



clei. Each object has been imaged in three continuum filters ('UV, 'blue' and 



'yellow') placed between the strong emission lines. The spectacular, high quality 
q ! images reveal a wealth of details about the structure of the host galaxies and 

their environments. Six of the nine galaxies in the sample are ellipticals with de 
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Vaucouleurs light profiles, one object has a well-defined disk component and the 
remaining two have marginal disks. Stellar populations of type II quasar hosts 
are more luminous (by a median of 0.3—0.7 mag, depending on the wavelength) 
and bluer (by about 0.4 mag) than are M* galaxies at the same redshift. When 
smooth fits to stellar light are subtracted from the images, we find both positive 
and negative residuals that become more prominent toward shorter wavelengths. 
We argue that the negative residuals are due to kpc-scale dust obscuration, while 
most positive residuals are due to the light from the nucleus scattered off interstel- 
lar material in the host galaxy. Scattered light makes a significant contribution 
to the broad band continuum emission and can be the dominant component of 
the extended emission in the UV in extreme cases. 

Subject headings: galaxies: active - - galaxies: structure - - polarization - 
quasars: general 

1. Introduction 

One of the key long-term questions in quasar studies is to identify just which properties 
of the galactic host (overall morphology? presence of bars? close neighbors? star-formation 
history? gas content?) are associated with activity in its nucleus. Observations of quasar 
hosts are made difficult by the small angular scales involved (1"~ 4.5 kpc for z = 0.3). 
Even for quasars with modest redshifts (z < 0.8), optical imaging from the ground requires 
excellent seeing conditions or the use of adaptive optics (Ronnback et al. 1996; Percival et 
al. 2001; Jahnke & Wisotzki 2003; Orndahl et al. 2003; Jahnke et al. 2004; see Bahcall et 
al. 1997 for earlier references). Imaging of quasar hosts can therefore greatly benefit from 
the superior angular resolution in space and was one of the major science drivers of the 
Hubble Space Telescope (HST; Bahcall et al. 1997; Boyce et al. 1998, 1999; McLure et al. 
1999; Hamilton et al. 2002; Dunlop et al. 2003; Pagani et al. 2003; Floyd et al. 2004). An 
additional difficulty present in ground-based and space-based observations is posed by the 
bright source in the center of the host galaxy - the quasar itself. In many cases the light 
from the quasar completely overwhelms the host galaxy, preventing its detection altogether 
(e.g., Hamilton et al. 2002). Indeed, the original definition of a quasar required a point-like 
morphology. 

In the context of lower luminosity Active Galactic Nuclei (Seyfert galaxies), much of 
the work on host galaxies has focused on "type 2" objects, Seyfert nuclei deeply obscured 
along our line of sight by matter close to the nucleus (Heckman et al. 1997; Malkan et 
al. 1998; Martini & Pogge 1999; Storchi-Bergmann et al. 2001; Gonzalez-Delgado et al. 
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2002; Kauffmann et al. 2003). There is strong evidence that type 2 Seyfert galaxies would 
be classified as ordinary, "type 1" objects if viewed from a different direction, where the 
obscuration is negligible (for reviews, see Antonucci 1993; Krolik 1999). The hosts of type 2 
Seyferts should then form a representative sample for the hosts of all Seyfert galaxies. The 
advantage of using type 2 Seyferts is that images of their hosts may be studied unencumbered 
by bright nuclei. Subtracting the light from the active nucleus is especially problematic at 
quasar luminosities, and a significant advantage could be gained by studying the hosts of 
obscured quasars. Investigation of obscured quasar hosts is the subject of the present paper. 

Although the direct line of sight to the nucleus is deeply obscured, nuclear light can 
contribute to the extended emission by way of scattering in the host galaxy (Antonucci & 
Miller 1985). The scattering efficiency (the ratio of the apparent luminosity due to scattering 
to the intrinsic luminosity) is very poorly known and certainly varies by a large factor from 
object to object, but some observations suggest that it can be as large as a few per cent 
(Zakamska et al. 2005). If the nuclear luminosity is ~ 100 x the host luminosity, a scattering 
efficiency of 1% doubles the apparent luminosity of the host. While such contrast between 
the nuclear and the host luminosities is rare in the V band (Bahcall et al. 1997), it might be 
common in the blue and near-UV, so we expect a significant contribution of the scattered 
light in these bands. Aided by polarimetric data and detailed morphological information, we 
attempt to separate the scattered light contribution from the stellar component of the host 
galaxy, but there are often severe limits on how well this can be done. 

In this paper we present results of the imaging of nine type II (obscured) quasars using 
the Advanced Camera for Surveys (ACS) aboard the HST, in combination with polarimet- 
ric data from ground-based observations. In Section 2 we describe the sample selection, 
observations and data reduction. We describe our tools for image analysis in Section 3. 
Individual objects are described in Section 4. We discuss our results in Section 5 and con- 
clude in Section 6. Throughout this paper, we use a cosmology with h = 0.7, Q m = 0.3 
and f^A = 0.7. Objects are identified as SDSS Jhhmmss.ss+ddmmss.s in Table 1 (e.g., 
SDSS J103951. 49+643004. 2) and as SDSS Jhhmm+ddmm elsewhere. Spectral features are 
identified using air wavelengths in angstroms (e.g., [OIII]5007). 

2. Sample selection and observations 

2.1. Sample selection 

The selection and optical properties of the parent sample of 291 type II Active Galactic 
Nuclei (AGNs) candidates were described by Zakamska et al. (2003), hereafter Paper I. 
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Briefly, these objects were identified in the spectroscopic database of the Sloan Digital Sky 
Survey (SDSS; York et al. 2000; Stoughton et al. 2002; Abazajian et al. 2004) as sources 
at redshifts 0.3 < z < 0.8 showing narrow emission lines with high-ionization line ratios 
characteristic of non-stellar ionizing continua. We used the [OIII]5007 emission line as a 
proxy for the intrinsic luminosity and estimated that those objects with [OIII]5007 emission 
line luminosities in excess of 3 x 10 8 L Q have intrinsic luminosities Mb < —23. This luminous 
subsample consists of about 150 objects which we classified as type II quasars (rather than 
Seyfert 2 galaxies). High bolometric luminosities (well in excess of the conventional quasar- 
defining value of 10 45 erg sec -1 ) of many of these objects and high densities of neutral gas 
along the line of sight were confirmed by infrared and X-ray observations (Zakamska et al. 
2004; Vignali et al. 2004; Ptak et al. 2006). 

At lower redshifts (0 < z < 0.33), type II AGNs from the SDSS have been studied 
extensively by Kauffmann et al. (2003) and Hao et al. (2005a,b). These samples were selected 
based on the emission line properties from the complete sample of SDSS galaxies (Strauss 
et al. 2002), but only a small number of objects can be classified as type II quasars based on 
the aforementioned [OIII]5007 luminosity criterion, with the majority of the objects being 
Seyfert 2 galaxies. 

We then restricted ourselves only to those of the type II AGNs from the samples by 
Zakamska et al. (2003) (Paper I) and by Hao et al. (2005a) which have L([OIII]5007)> 1O 9 L 
(i.e., well above the quasar luminosity criterion) and which were radio-quiet (Zakamska et 
al. 2004). Since the equivalent widths of the emission lines in type II quasars can be very 
high (above ~ 1000A, Paper I), these lines can make a significant contribution to the broad- 
band flux. To avoid contamination from the strongest lines ([Oil] 3727, [OIII]4959,5007, H/3, 
Ha) to the HST images, we selected out objects from two narrow redshift ranges so that 
three ACS broad-band filters could be placed between these lines. We selected three objects 
at z ~ 0.25 from the sample by Hao et al. (2005a) and six objects at z ~ 0.4 from the 
sample in Paper I. SDSS spectra of these objects are shown in Figure 1 together with the 
wavelength coverage of the filters that we used in our ACS observations. Using the [OIII]5007 
line luminosities, we estimate that the intrinsic luminosities of the objects described in this 
paper are about —24 > M B > —26. 

Throughout this paper, we refer to the filters used in our observations as 'yellow', 'blue' 
and 'UV, in order of decreasing wavelength. The 'yellow' band has an effective rest-frame 
wavelength of 5700— 5800A for all objects in our sample, i.e., somewhat redder than but 
close to the rest-frame V band, and the 'blue' and 'UV bands are longward and shortward 
of [Oil] 3727, respectively. 
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2.2. HST observations and data reduction 

All observations pertaining to our program were obtained in the period July 2003 - May 
2004. Each object was imaged with the Wide Field Channel (image scale 0.049" pixel -1 ) of 
the ACS in three filters, as summarized in Table 1 and illustrated in Figure 1. For each filter, 
four exposures were taken, each lasting 540-610 sec, in a four-point dither pattern, or overall 
observing time of one orbit per object per filter. Data reduction consisted of rejecting cosmic 
rays from each of the four flat-field-corrected exposures, applying a geometric distortion 
correction and then combining the four images. All data reduction was performed using the 
MultiDrizzle routine (Koekemoer et al. 2002) available as part of the HST archive pipeline. 

Using the standard ACS photometric calibration described in the ACS data handbook 
(2005), we converted the units of the reduced data (counts/sec) into spectral flux densities F v 
and AB magnitudes itlab = —2.5 logF u — 48.60, where F v is in units of erg sec -1 cm -2 Hz -1 . 
The absolute AB magnitudes were calculated at rest-frame effective wavelengths X bs/ (1 + z ) 
as Mab = m AB — 5 log(Di,/10 pc) + 2.5 log(l + z), where is the luminosity distance and 
\ bs are the effective wavelengths in the observer's frame. The values \ b s / (1 + z) are listed 
in Table 1. 

The images taken in three bands were rebinned to the common grid using point sources 
in the field and combined to produce a color-composite image using the asinh stretching 
and color-combining routine by Lupton et al. (2004). The resulting images are presented 
in Figures 2-3. We used identical image combining parameters for the sources with similar 
redshifts that used similar or identical filters (one group of five objects and another group of 
three objects, Table 1). As a result, images from the same group can be directly compared 
with one another. For example, the yellow color of SDSS J0139+6430 compared to the red 
color of SDSS J0123+0044 reflects a difference in their stellar populations and is not an 
image processing artifact or the effects of K-corrections. The images in Figure 3 (group 
2 objects) have a rather different color scheme from those in Figure 2. This difference is 
due to using the ramp filter FR459M as the UV filter for the group 2 objects; this filter is 
significantly narrower than all other filters in our program (Figure 1), resulting in about a 
factor of 10 reduction in sensitivity. Finally, SDSS J1323— 0159 was not grouped with any 
other sources because it is the only object that uses F775W as its yellow filter. 

2.3. Optical polarimetry and spectropolarimetry 

In Table 2 we summarize all polarimetric and spectropolarimetric data available for the 
objects observed with HST. Out of nine objects, broad-band polarimetry is available for two 
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objects and spectropolarimetric data are available for five objects. Of the latter, three objects 
(SDSS J1039+6430, SDSS J1323-0159, and SDSS J1413-0142) were previously reported by 
Zakamska et al. (2005) (hereafter Paper II), and the results of these observations are given 
in Table 2. In this section, we describe spectropolarimetric and polarimetric observations of 
four more objects. 

Spectropolarimetry was obtained for two objects, SDSS J0920+4531 and SDSS Jl 106+0357, 
on December 19, 2004. We used the CCD Spectropolarimeter (SPOL, Schmidt et al. 1992b) 
at the 6.5m Multiple Mirror Telescope (MMT). All observations used a low-resolution grat- 
ing providing spectral coverage of 4100-8200A. With seeing of more than 2" (full width at 
half maximum), an entrance slit of 1.5" width was used, resulting in a spectral resolution of 
~22A. Each object was observed for 6400 sec. The absolute flux calibration was performed 
using the original SDSS spectra (whose spectrophotometric calibration is accurate to 5%, 
Abazajian et al. 2004) rather than photometric standards because seeing was highly variable 
on the night of the observations. The instrumental polarization was verified during the runs 
to be < 0.1% through observations of unpolarized standard stars taken from Schmidt et 
al. (1992a), and the polarization due to intervening Galactic dust is < 0.2% and < 0.5% 
for SDSS J0920+4531 and SDSS J1106+0357, respectively (estimated using the Galactic 
reddening values by Schlegel et al. 1998 and the maximum interstellar polarization per unit 
reddening by Serkowski et al. 1975). Other details of instrument configuration, calibration 
and data reduction are the same as described in Paper II. 

The results of the spectropolarimetric measurements are presented in Figure 4. The 
degree of polarization is quantified by the rotated Stokes parameter g(A) = [<5acos26* + 
U\ sin 29}/ F\, where 9 is the mean polarization angle of the source in the 5000-8000A band- 
pass, F x is the total optical flux density, and Qx and Ux are Stokes parameters. Unlike 
the conventional polarization P(A) = y/Q\ + U^/Fx, the value q(X) is not positive def- 
inite and has a normal error distribution. The polarization position angle is defined as 
tan 20(A) = Ux/Q\ and is measured in degrees East of North. In Figure 4, the degree of 
polarization, the polarization angle and the polarized flux q(X) x Fx have been heavily binned 
in wavelength because of the low signal-to-noise ratio at the original spectral resolution. In 
Table 2, we summarize spectropolarimetric data in the form of broad-band measurements 
('UV, 'blue' and 'yellow') obtained by error- weighting spectral data within the corresponding 
ACS filters' bandpasses. 

Broad-band polarimetry was obtained for 11 objects from Paper I in September 2004 
(see Appendix A for the complete dataset). Among them were two objects from our HST 
program (SDSS J0123+0044 and SDSS J2358-0009). Observations were conducted using 
SPOL in its imaging mode at the Bok 2.3-m reflector on Kitt Peak, AZ. The instrumental 
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configuration and data processing were similar to those discussed by Smith et al. (2002). In 
brief, a semiachromatic wave plane is placed in the beam and rotated through 16 positions 
to determine the linear Stokes parameters, Q and U, using exposure times for each position 
of the plate of 200-240 sec, with a total integration time per object of 3200-3840 sec. All 
observations were through a V filter (4800— 6000A) to prevent [OIII]5007 from being included 
in the bandpass. The resulting broad-band polarization values (again quoted as rotated 
Stokes parameters), position angles and their uncertainties are listed in Table 2. In both 
objects polarization is detected at less than 2a confidence level. 

3. Image analysis and morphology 

Our morphological analysis of the HST data is aimed at determining (i) whether these 
galaxies can be described as disks, disks with bulges, or pure ellipticals, (ii) whether there 
is any evidence for a non-stellar contribution to the HST images, (iii) whether there are 
significant morphological distortions of the hosts indicative of recent or ongoing interactions, 
and (iv) the luminosities and colors of the stellar populations of the host galaxies. We start 
by performing aperture photometry of our objects in Section 3.1 and then proceed to derive 
model photometry. We summarize the definition of the Sersic profile in Section 3.2. In 
Section 3.3 we describe one-dimensional (1-D) surface brightness profiles of the galaxies in 
our sample. We then use the results of the 1-D analysis in fitting smooth two-dimensional 
(2-D) elliptical models in Section 3.4. 

3.1. Aperture photometry 

For a given object, we used elliptical apertures with ellipticities, position angles and 
semi-major axes determined by the outer isophotes in the yellow image. The same aperture 
was used in all three bands for this object. We tried several aperture sizes and chose the 
ones which were large enough so that further increase in aperture size did not increase the 
total enclosed flux by more than 5%, while at the same time avoiding close companions. The 
aperture sizes were in most cases somewhat larger than the outermost detected isophote (our 
typical surface brightness sensitivity is 25—26 mag/arcsec 2 in the yellow band). We estimate 
that the resulting aperture magnitudes (Table 3) are accurate to better than 0.1 mag. 
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3.2. Sersic brightness profile 

If the isophotes of a galaxy are elliptical with the same ellipticity e = 1 — (semi-minor 
axis/semi- major axis), the same position angle on the sky and the same center, its photom- 
etry can be completely described by the 1-D surface brightness profile I(r), where r is the 
semi-major axis of the isophote. The Sersic (1968) profile is a three-parameter function: 

log J(r) = a - air l/n % (1) 

where n s is referred to as the Sersic index. The surface brightness profiles of many galaxies 
can be described by the Sersic function or a superposition of two functions with different n s 
(Binney & Merrifield 1998). Elliptical galaxies and bulges of spirals typically have n s ~ 4.0 
(de Vaucouleurs 1948), while disk components typically have an exponential profile, i.e., 
n s ~ 1.0 (Patterson 1940; Freeman 1970). In equation (1), parameter a\ determines the scale 
length of the galaxy. For the exponential and the de Vaucouleurs profiles this parameter 
can be related to the conventionally used exponential scale length and half-light radius, 
correspondingly: Rd = 0.434/ai and R e = (3.33/ai) 4 . 

Assuming that the profile continues to infinity, the total flux emitted by a galaxy with 
a 1-D surface brightness profile given by equation (1) and with ellipticity e is 

F={l-e)2vW \ ns ^ 2n j , (2) 
V ; (ailnlO) 2 "* V ; 

If n s is an integer (as it is for the de Vaucouleurs and the exponential profiles), then T(2n s ) = 
(2n s - 1)! 



3.3. 1-D brightness profiles 

For each object and each band, we use the standard IRAF package ellipse to fit isophotes 
with ellipses, without regard to whether the light is due to the stellar or the scattered 
component. Isophotes are fitted at pre-defined semi-major axes using the iterative algorithm 
described by Jedrzejewski (1987). For each semi-major axis, the output includes the best-fit 
surface brightness, center, ellipticity, position angle and other parameters. We then plot the 
surface brightness as a function of the semi-major axis, without taking into account changes 
in the centers, ellipticities or position angles from one isophote to another. The resulting 
1-D profiles are shown for all objects in all filters in Figure 5 both as a function of r and r 1 / 4 . 
In most cases (especially in the yellow and blue bands) the models constructed based on 
best-fit elliptical isophotes reproduce the observed brightness profile very well (better than 
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1%), but in some cases (typically in the UV) significant residuals are present; examples are 
shown in Figure 6. 

The point-spread function (PSF) suppresses the central flux and redistributes it outward, 
causing the flattening of the inner parts of each profile (especially apparent in the r 1 / 4 — 
log J(r) plots in Figure 5). The PSF also circularizes the isophotes and slightly increases the 
effective radii. We explicitly incorporate the PSF in our 2-D fits (Section 3.4). 

Although the constructed 1-D brightness profiles do not take into account the changes 
in the ellipticities, position angles or positions of the centers of the isophotes with semi- 
major axis (and neither do they discriminate between scattered and stellar components), 
they provide some insights into the structures of the galaxies and a starting point for further 
analysis (Section 3.4). Figure 5 shows that the de Vaucouleurs profile fits many of the galaxies 
very well. Strong evidence for an exponential component is present only in SDSS J1106+0357 
(which indeed looks like a spiral galaxy, Figure 3). The 1-D profiles of SDSS J0920+4531 
and SDSS J1413— 0142 show changes in slope, possibly indicating a transition from a bulge 
to a disk component (although these changes do not occur consistently in all three bands). 
In some objects, e.g. SDSS J1323— 0159, there is a clear color gradient that can be seen as 
different slopes of the brightness profiles in three bands - in this case, the UV light is more 
centrally concentrated than blue or yellow. 

A number of other features can be identified in the 1-D brightness profiles. For example, 
the prominence of wiggles, such as those seen in the 1-D profiles of SDSS J0920+4531 and 
SDSS J1413— 0142, increases toward the blue, suggesting that they are due to scattered light 
or patches of star formation. Some centers are clearly underluminous compared to the best- 
fit smooth profiles (e.g., SDSS J1323— 0159), while others show a strong central excess of 
light despite the smoothing due to the PSF (e.g., SDSS J1301— 0058). We will argue that flux 
deficiencies (seen as negative residuals when the best-fit 2-D model is subtracted) are caused 
by dust obscuration, whereas the interpretation of flux excesses (seen as positive residuals) 
may be due to scattered light, patches of star formation or a superposed companion galaxy. 

3.4. 2-D brightness profiles 

In Paper II, we demonstrated that at least in some cases light emitted by the central 
AGN and scattered off interstellar material in the host galaxy contributes a significant portion 
of the total flux and therefore should be present in the HST images. In Figures 2-3 we show 
the orientations of the brightest UV spots or cones and the orientation of the polarization 
vector (where polarimetric data are available). In some cases (such as SDSS J1323— 0159 and 
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SDSS J1413— 0142 described in Paper II) the scattering regions can be identified based on 
their morphological similarities to well-studied nearby examples (Pogge & de Robertis 1993), 
and this identification agrees with the orientation predicted by spectropolarimetric measure- 
ments. In most cases, the polarization position angle is orthogonal to the line between the 
centers of the galaxies and the presumed scattering regions, but there are exceptions (such 
as SDSS J0920+4531, where the brightest UV spot does not seem to be responsible for 
the polarization) and complications (such as SDSS J0123+0044, where numerous spots are 
detected). In addition, the polarization signal is suppressed when scattering occurs along 
multiple directions, partially canceling the polarization vectors; in this case the scattering 
regions cannot be unambiguously identified from the polarimetric data. 

The major complication of our analysis is distinguishing between scattered and stellar 
components. In the absence of detailed imaging polarimetry some assumptions must be 
made about the two components to allow their separation. We fit our galaxies with a 2-D de 
Vaucouleurs model (elliptical galaxy), or a sum of a de Vaucouleurs model and an exponential 
model (disk galaxy with a bulge) in the cases where the presence of a disk is suggested by 
the form of the 1-D profile (Section 3.3). We assume that the stellar component of the host 
galaxy is well represented by these simple models, whereas the residuals are comprised of 
scattered light, isolated regions of star formation and dust obscuration. We comment on the 
robustness of the scattered light identification for each object in Section 4. For each object, 
the images in the three bands are modeled independently. With the exception of UV-band 
images of group 2 objects, the galaxies are typically detected well outside two effective radii, 
and out to five or so effective radii in some cases, allowing robust fitting of the outer parts. 
Our adopted models are certainly overly simplified, but the quality of our data and the 
presence of strong irregular morphological features prevents the use of more complex models 
with a larger number of free parameters (for an extensive library of galaxy models see, e.g., 
Peng et al. 2002). 

We convolve all models with the PSF before fitting the data. We used TinyTim 1 to 
generate PSFs for all filters in our observations. For the broader filters, the PSF has some 
dependence on the assumed spectral energy distribution within the filter. In these cases, we 
used an approximation to the spectrum in the form F\ =const. The size of the PSF kernel 
used in the convolution is 3". 

In the simplest case of a de Vaucouleurs model, there are six fitting parameters: total 
normalization ao, slope ai, ellipticity, the position of the center (two coordinates), and the 
position angle of the major axis. We assume concentric, aligned isophotes with ellipticities 



1 http://www.stsci.edu/software/tinytim 
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that are constant with the distance from the center. In the bulge+disk model, there are 12 
parameters (6 for each of the Sersic components) that are allowed to vary independently. 
The mean value of the sky is kept constant during 2-D fitting after it is determined by fitting 
a Gaussian to the distribution of pixel values for a patch of blank sky in the same image. 
The value being minimized to find the best fit is Ylii^i ~ Mi) 2 /<jf, where summation is over 
all unmasked pixels, is the flux of the i-th pixel, Mj is the model value of the same pixel 
and a, is an estimate of the flux error per pixel that combines the Poisson noise and the 
read-out noise. 

The most challenging part in fitting the 2-D profiles to the data is to identify and mask 
out regions with significant deviations from the Sersic profiles. Polarization data help identify 
the major scattering regions as demonstrated in Figures 2-3. We used an iterative procedure 
in which the fitting residuals are used to construct a mask for the next round of fitting, and 
the procedure is repeated until the residuals very closely resemble the mask. In some cases 
the results of the fitting are sensitive to the assumed mask, leading to large uncertainties 
in structural parameters such as half-light radii or ellipticities. However, the errors in these 
parameters are highly covariant, with the result that the uncertainty in the measurement of 
the integrated flux is rather small. We estimate errors by running a series of 2-D fits with 
different masks, selecting acceptable fits and looking at the range of structural parameters 
and luminosities that they yield. 

In Table 4 we list the best-fit structural parameters that resulted from our 2-D modeling. 
The absolute magnitudes listed in Table 4 are obtained using equation 2 ('model magnitudes' 
hereafter), i.e., assuming that the best-fit profiles extend to infinity. 

Whenever our fitting procedure results in negative residuals, these residuals typically 
increase toward shorter wavelengths, so it is natural to attribute them to dust obscuration. 
In these cases, we calculate the ratio of observed to model flux in the yellow band (which is 
close to the rest-frame V band in all cases) and correct it to the V band using an extinction 
law in the form A\ oc 1/A (e.g., Cardelli et al. 1989), so that all values of extinction are 
given corrected to the rest-frame K-band. These corrections are much smaller than the 
uncertainties in the model. 

In order to compute the ratio of the scattered to stellar light listed in Table 4, we first 
calculate the total stellar luminosity using equation 2 and then perform aperture photom- 
etry on the positive residuals that we confidently attribute to scattered light based on the 
polarization measurements. In Section 4, we explicitly describe which regions were included 
for each object. Our procedure is likely to underestimate the fractional importance of the 
scattered light because we include only those positive residuals that are unambiguously due 
to scattering, and because low surface brightness scattered emission may disappear in the 
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fit, may be underestimated because of dust obscuration or may simply remain undetected. 

4. Individual objects 

4.1. SDSS J0123+0044, z = 0.399 

This object is represented by a single de Vaucouleurs profile quite well, as shown in 
Figure 7. Residuals become stronger toward shorter wavelengths, and some dust obscuration 
(Ay < 0.2 mag) is present near the center. The nature of the residuals is ambiguous. Some 
of the bright UV spots in the composite image (Figure 2; also seen as residuals in Figure 
7) may be scattering regions, since the line connecting the brightest spot with the nominal 
center of the galaxy is orthogonal to the polarization position angle. This interpretation 
should be regarded as somewhat uncertain, since the polarization was measured at less than 
2o- confidence level. On a larger scale, the galaxy shows a faint tidal tail-like structure 
(Figure 8) indicating possible interactions. Therefore, it is likely that at least some of the 
residuals are star forming regions, especially those to the east of the nucleus that connect to 
the tidal tail. In our estimate of the scattered light contribution given in Table 4, we do not 
include any of the positive residuals to the east of the nucleus. 

4.2. SDSS J0920+4531, z = 0.402 

This object shows some evidence for an exponential component in its 1-D profile, so we 
fit it with a bulge+disk combination, producing rather poor fits spanning a large range of 
possible fitting parameters depending on the chosen masks. The presence of dust obscuration 
and strong positive residuals with complicated morphology make our best fits particularly 
uncertain in this object. Dust obscuration is estimated at Ay < 0.5 mag; it can be seen as a 
red spiral-like structure in the color-composite image in Figure 2 to the east of the nucleus, 
suggesting that the eastern side is closer to the observer. We estimate that the uncertainty 
in the measurement of the total stellar magnitude in this object is about 0.3 mag, compared 
to <0.1 mag in the objects with better fits. This is the only object in the sample in which 
the brightest UV region in the color-composite image (the one to the north-west of the 
nucleus schematically shown as a grey cone in Figure 2) is not orthogonal to the measured 
polarization position angle. In fact, the presence of at least three morphologically irregular 
companions within 25 projected kpc and the morphology of the UV region suggests that it 
might be a superposed companion galaxy with on-going star formation. We therefore do not 
make an estimate of the scattered-to-stellar ratio for SDSS J0920+4531, since the scattering 
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regions cannot be confidently identified. 

4.3. SDSS J1039+6430, z = 0.402 

We model the stellar light in this object with a single de Vaucouleurs profile, as shown 
in Figure 9. We interpret all of the positive residuals as scattered light, taking into account 
the very high measured polarization in this source, which reaches 17% in the UV (Paper II). 
Some dust obscuration (A v < 0.3 mag) may be present to the north of the nucleus. After 
computing the scattered-to-stellar ratio using aperture photometry of the positive residuals, 
we correct the observed values of polarization in all three bands (Table 2) for the contri- 
bution of the stellar light (assumed unpolarized) . We find the polarization of the scattered 
component to be 26±3%, independent of wavelength within the errors. An additional cor- 
rection can be made to account for the finite width of the scattering cone: particles at 
different positions within the cone produce polarization at somewhat different polarization 
position angles, leading to partial cancellation of polarization. If the opening angle of a 
homogeneously-filled scattering cone is 9q, the effect reduces the intrinsic polarization by a 
factor of sm9 /9 . For SDSS J1039+5430, 6 ~ 50°, so our final estimate for the intrinsic 
polarization is 30%. 

The spectacular large-scale (10 kpc in projection) scattering region to the south of the 
center was discussed in Paper II. Based on the large size of the scatterer, we argued that 
dust rather than electron scattering produces the high observed polarization in this object. 
This argument is further supported by the one-sidedness of the scattering region; although 
its counterpart to the north of the nucleus is detected and can be seen as a bluish spot on 
the color-composite image (Figure 2) and in the residuals (Figure 9), it is much fainter. The 
difference in the brightness of the northern and the southern spots can be explained due 
to the preferentially forward scattering of dust particles, if the axis of the scattering cones 
points somewhat toward the observer in the southern part. 

4.4. SDSS J1106+0357, z = 0.242 

This object is the largest and best-measured bulge+disk galaxy in our sample. In the 
yellow and blue bands, significant deviations from the best fitting bulge+disk are present 
(Figure 10). Obscuration is present both in the outer disk (Ay < 1 mag) and in the central 
parts of the galaxy (A v < 0.6 mag). In the outer disk the obscuration is preferentially 
confined to the northern (image top left) side, suggesting that this side is closer to the 
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observer. The scattering region is well defined and has a classical biconical shape, and it is 
accurately orthogonal to the measured polarization position angle. Although the scattering 
regions are detected out to 1.5 kpc from the center, they are dwarfed by the galaxy which 
has a disk exponential scale well above 10 kpc (Table 4). The disk component is not detected 
in the UV-band image, presumably because of the low sensitivity of the ramp filter used in 
this observation, and in Table 4 we give a lower limit on the UV-band stellar luminosity of 
the bulge component only. In the yellow image, at least three lumps can be seen outside the 
main disk but within 25 kpc (projected) of the center of this galaxy. A faint arm (seen in 
Figure 10) connects the outskirts of the disk with an extended companion to the north-east 
(image left) of the center. If this feature lies in the plane of the galaxy and if the observed 
ellipticity of about 0.4 is entirely due to inclination, then this arm follows a nearly circular 
arc around the center of this galaxy, with a radius of 25 kpc. 



4.5. SDSS J1243-0232, z = 0.281 

This object is very well represented by a de Vaucouleurs profile. However, there is a 
significant central flux excess over the best-fit de Vaucouleurs component within the few 
central pixels (Figures 5 and 11). The excess is very blue in color, approximately Fx oc A~ L5 , 
but whether this excess is due to a circumnuclear starburst or to sub-kpc scattering is not 
clear. This excess is excluded from the calculation of the model magnitudes. 



4.6. SDSS J1301-0058, z = 0.246 

We start by analyzing the largest component in this guitar-shaped interacting galaxy. 
As in SDSS J1243— 0232, the 1-D profile of this galaxy shows a clear central excess over the 
de Vaucouleurs fit to the outer parts of this object. A combination of two de Vaucouleurs 
profiles (one representing the central excess and one representing the outer parts) provides the 
best fit to the 2-D images. The more extended de Vaucouleurs component is clearly of stellar 
origin, but the nature of the more compact component (R e ~ 0.6 kpc) is unclear. As for 
SDSS J1243— 0232, we only include the data for the more extended of the two components 
in the calculation of the luminosity of the stellar component (listed in Table 4). Some 
residuals are present when the best fit is subtracted. The orientation of these residuals is 
schematically shown in Figure 3, but since no polarization data are available for this object, 
we cannot determine whether if these residuals are caused by scattered light. The companion 
galaxy to the south-west (lower left) of the main object is well-represented by a disk+bulge 
combination and shows UV bow-like residuals in the center that could also be interpreted 
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as scattering cones. 

It is not clear which of the two objects hosts the type II quasar. The separation between 
the galaxies is 1.3", so that they both comfortably fit into the 3" diameter SDSS fiber 
aperture. Both galaxies have irregular central parts, so given the range of morphologies of 
the scattered light in our sample, either of the two objects (or both) could be the actual 
AGN host. 



4.7. SDSS J1323-0159, z = 0.350 

The outer parts of this galaxy are well-fit by a de Vaucouleurs profile (Figure 12), but 
there are significant non-axisymmetric negative and positive residuals in the center, making 
2-D fitting of this object quite challenging. These difficulties are manifested by the range of 
the effective radii allowed by our fits shown in Table 4. Nevertheless, the uncertainties in 
the measured absolute magnitudes are quite small (<0.1 mag) in the yellow and blue bands. 
In the UV, the uncertainties increase to 0.2 mag because the stellar light is more compact 
in this band and the fit depends more sensitively on the assumed central mask. 

This object was first presented in Paper II, where we argued, on the basis of morphology, 
that the UV X-shaped structure (seen as positive residuals in Figure 12) is a biconical 
scattering region. This interpretation is supported by the fact that the axis of this structure 
is orthogonal to the polarization position angle. The negative residuals can be seen as a red 
streak going vertically through the center of the color-composite image of this object (Figure 
2) and can be interpreted as a kpc scale dust lane with the largest observed extinction about 
Ay ~ 1.2 mag. 

4.8. SDSS J1413-0142, z = 0.380 

The 1-D surface brightness profiles of this galaxy show evidence for two breaks (at 
1.5 kpc and at 4.5 kpc in semi-major axis), and the profile between the breaks is roughly 
exponential. In our 2-D fitting procedure, we fit this object with a disk+bulge combination. 
As for SDSS J0920+4531, this results in rather poor fits. 

This object was first presented in Paper II and shows two bright conical structures 
(schematically represented by grey cones in Figure 2) roughly orthogonal to the polarization 
position angle, although not quite as symmetric around the center as the more classical 
cones of SDSS J1323— 0159. In addition to these prominent structures, excess light is seen 
in the residuals and in the color composite to the south and to the north of the nucleus. 
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Furthermore, residuals show that some extinction is present in the central region. As a result, 
the fits for the bulge are highly uncertain, especially in the UV band. In our two- component 
model of the galaxy, the disk is offset relative to the bulge by about 0.4 - 0.65 kpc from the 
center along the semi-major axis, roughly toward the south of the nucleus (this asymmetry 
can be seen in Figure 2). In the estimate of the scattered flux, we include all of the positive 
residuals because of their conical morphology. 

4.9. SDSS J2358-0009, z = 0.402 

The 1-D profile of this object suggests a de Vaucouleurs profile, which produces a good 
2-D fit, but with strong positive central residuals with close-to-zero ellipticity. 2-D modeling 
of the UV image produced a significantly larger effective radius than that suggested by the 1- 
D profile, and therefore the UV effective radius is not listed in Table 4. The ambiguity is due 
to the fact that most of the detected UV light is masked out from the fit, as it is part of the 
central excess rather than the extended stellar light. It is not clear whether these residuals 
are due to scattered light or a circumnuclear starburst. In addition, a narrow UV cone 
extending to the north-east of the nucleus and orthogonal to the polarization position angle 
(Figure 2) is well-detected in the residuals, and we interpret this feature as scattered light 
and include it in Table 4, although we caution that the polarization is detected at less than 
2c confidence level in this object. Curiously, the UV cone is aligned with a tidal structure 
connecting SDSS J2358— 0009 to the companion galaxy, so another possible interpretation 
is a stretched star-forming region. 

Because the UV feature appears jet-like, we also test the possibility that it represents 
synchrotron radiation from an energetic outflow, using 1.4 GHz radio data from the FIRST 
survey (Becker et al. 1995; White et al. 1997) and the NVSS (Condon et al. 1998). There 
is neither a FIRST nor an NVSS detection at the position of SDSS J2358-0009, and using 
the rms noise in the FIRST image of the field we place a 5a upper limit of 0.7 mJy on the 
flux from this source (assuming a point source). From our HST data, a lower limit on the 
UV flux density from the conical feature (excluding the central part where the morphology 
of the positive residuals is uncertain) is F\ ~ 4.2 x 10~ 19 erg sec -1 cm -2 A -1 . Therefore, the 
radio-to-optical index of this feature a (defined as F opt / F ra< n = {v pt/ ^radio) a ) is constrained 
to be > —0.6 by our data. While values around a = —0.5 are not uncommon for the total 
emission (including jets, lobes, radio cores, host galaxy) from radio-loud sources, the jets 
have spectra with typical values of a around —1 (Meisenheimer et al. 1996; Sambruna et al. 
2004; Jester et al. 2005). Therefore, the observational limit on a makes SDSS J2308-0009 
an unlikely jet candidate. 
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The galaxy is clearly experiencing some type of interaction, showing large-scale tidal 
tails (Figure 13) as well as smaller scale shell-like structures to the east of the nucleus. 
The interacting companion to the north-east (at the projected distance of 29 kpc) is well 
represented by a sum of disk and de Vaucouleurs component. 

5. Discussion 

5.1. Stellar populations of the type II quasar hosts 

In this section we discuss the luminosities and colors of the stellar population of the type 
II quasar hosts we have imaged using the model magnitudes and characteristic radii (Table 4) 
that resulted from our 2-D fitting procedure. To compare the luminosities of the hosts in our 
sample to those of the galaxies from the general population at the same redshift, we calculate 
the so-called M* magnitudes defined as the position of the break of the Schechter (1976) fit to 
the luminosity function of galaxies. The number density of galaxies with magnitudes brighter 
than M* decreases roughly exponentially with luminosity, so M* quantifies the luminous end 
of the luminosity function; Bahcall et al. (1997) estimate that the average luminosity of a 
field galaxy is 1.8 mag fainter than M*. We used the luminosity function of galaxies from 
Blanton et al. (2003) at five optical wavelengths. We corrected their values of M* to match 
the redshifts of our objects using their own best-fit redshift evolution parameters. We then 
used linear interpolation to calculate M* values at the effective rest-frame wavelengths of our 
observations listed in Table 1. The resulting M* values are presented in Table 4 in parentheses 
next to the model magnitudes of our galaxies. Although the M* values by Blanton et al. 
(2003) have nominal accuracies of < 3%, our calculation of M* involves extrapolation of 
their results to higher redshifts (z < 0.4) than those of the galaxies used in constructing the 
luminosity functions (z < 0.22), so we quote only three significant digits for our computed 
M* in Table 4. 

In Figure 14a we compare the absolute model magnitudes of the stellar components of 
out type II quasar hosts with M* values. This comparison shows that in the yellow band 
five of the nine objects are within 0.3 mag of M*, while the remaining four are 0.7—1.4 mag 
more luminous. The median difference is (Mh ost — M*) mei i = —0.3 mag in this band. The 
luminosity difference between type II quasar hosts and field galaxies is more pronounced 
in the blue band, in which type II quasar hosts are up to 2 mag more luminous than the 
M* values, with a median difference of 0.7 mag. Similarly, the hosts of unobscured quasars 
are often found to be significantly more luminous than are field galaxies, but quantitative 
comparison is complicated. For example, Bahcall et al. (1997) find (M host — M^) med = mag 
and Hamilton et al. (2002) find a —0.65 mag difference (both these values were obtained 



-18- 



in the V band which is close enough to our yellow band to allow a direct comparison, 
and both values include only radio-quiet quasars). The differences in the host luminosities 
computed by different groups are at least in part due to the modeling difficulties associated 
with subtracting the bright central source, but may also be related to the differences in 
the luminosities of quasars in different samples, if host and quasar luminosities are strongly 
correlated. 

In Figure 15 we show type II quasar nuclear and host luminosities, as compared with the 
same values for unobscured quasars from the literature. Because the quasars in our study 
are obscured in the optical, their intrinsic My values are computed using the [OIII]5007 
line luminosity (Paper I); the error bars reflect the la uncertainty due to the scatter in 
the [OIII]5007— My correlation for unobscured quasars. Our yellow band is close to the 
V-band, and the K-corrections required to obtain V-band host magnitudes from the yellow- 
band magnitudes listed in Table 4 are less than 0.1 mag. For unobscured quasars, the PSF 
of the bright central source has to be included in the 2-D modeling, and the difference 
between host/nuclear luminosities obtained by Bahcall et al. (1997) and Hamilton et al. 
(2002) for identical objects (solid lines in Figure 15) can be taken as a measure of the 
2-D modeling uncertainties. Any 2-D fitting procedure based on x 2 minimization will be 
somewhat degenerate in separating the nuclear flux from the host flux, but will keep the 
total flux roughly constant. Indeed, most of the solid lines follow the contours of constant 
flux (examples are shown in dotted lines) rather closely. 

From Figure 15, we find that type II quasar hosts occupy a similar range of luminosities 
to the hosts of luminous quasars. We also find that the uncertainties in the positions of type 
II quasars on the nucleus/host diagram are comparable to those of unobscured quasars. The 
large uncertainties in the nuclear luminosities of type II quasars are counter-balanced by our 
~ 0.1 mag photometry of the hosts. 

Figure 14b shows the distribution of the colors of type II quasar hosts relative to those of 
the M* galaxies. Type II quasar hosts have UV— yellow colors similar to those of M* galaxies, 
but they are significantly bluer in the blue— yellow color, with a median color (relative to 
M* values) of —0.4 mag, reaching —0.8 mag for SDSS J1039+6430. There is no evidence 
that the colors of the disk+bulge systems in our sample are different from the colors of the 
ellipticals. 

The values of the colors can be used to place constraints on the ages of the stellar 
population. For example, for the six objects at redshift around 0.4 (i.e., excluding group 2 
objects) the colors can be directly compared with one another as the effective wavelengths 
of the observations are very similar in all six cases; the median colors for these objects are 
0.8 mag (blue— yellow) and 1.8 mag (UV— yellow). To estimate the age corresponding to 
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these colors, we use the stellar population synthesis models of Bruzual & Chariot (2003) of 
a passively evolving instantaneous starburst and produce a series of spectra as a function 
of age. We then calculate the model colors using these model spectra and the effective 
wavelengths of our observations listed in Table 1. For a fixed age, of all the parameters 
probed by Bruzual & Chariot (2003), such as different evolution prescriptions, different 
libraries of UV stellar spectra and different initial mass functions, the one that affects the 
colors the most is metallicity. If the initial metallicity is assumed to be solar, then the median 
colors of our galaxies indicate an age between 3xl0 8 and 10 9 years. If the initial metallicity 
is 20% solar, then the inferred age is around 10 9 years. For SDSS J1039+6430, the bluest 
object in the sample, the age estimates lie between 10 8 and 3 x 10 8 years. Our age estimates 
were obtained under the assumption that the entire galaxy underwent an episode of star 
formation, so if only a fraction of the mass of the galaxy was involved, then our estimated 
ages for the starburst should be regarded as upper limits. 

Our estimated ages suggest that A-type stars are an important contribution to the type 
II quasar hosts. They are best probed by our blue band observations which cover the spectral 
region near the Balmer break. In Paper II we demonstrated, using the high signal-to-noise 
ratio observations of the continuum emission, that the stellar light of SDSS J1039+6430 is 
dominated by an A-star population. 

The effective radii listed in Table 4 (exponential scales in the case of the SDSS Jl 106+0357 
disk component) can give us a handle on the color gradients in the stellar populations. The 
accuracy of the measurement of the a± parameter is much worse than the accuracy of the 
model photometry (Section 3.4), and furthermore the effective radius is oc a^ 4 for a de Vau- 
couleurs profile, so we estimate that the accuracy of the measurement of the effective radii is 
no better than 20%, limiting us to a qualitative description. Focussing only on the extended 
emission that is well-modeled by the de Vaucouleurs profile (or exponential in the case of 
SDSS Jl 106+0357) and disregarding the central excesses of ambiguous origin, we notice that 
our sample presents a mix of galaxies that become bluer outward (e.g., SDSS J1243— 0232) 
and those with the color gradient of the opposite sign (e.g., SDSS J1323— 0159). 



5.2. Scattered light 

Free electrons and dust particles in the interstellar medium of the host galaxy can scatter 
nuclear light. Distinguishing between these two possibilities can be rather difficult, as has 
been demonstrated by many authors (see Paper II for a detailed description and references), 
and it seems that both mechanisms operate in practice (e.g., Kishimoto 1999; Hines et al. 
2001). For our sample, we find strong evidence for dust scattering in SDSS J1039+6430 
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(Papcr II and Section 4.3), but similar arguments are not conclusive when applied to other 
objects. Both for electron scattering and dust scattering the number density of particles 
required to account for the observed amount of scattered light can be computed from the 
geometric dimensions of the scattering regions, supplemented by a guess about scattering 
efficiency (Paper II). Interestingly, the required number density turns out to be n e > 1 cm~ 3 
for electron scattering or n H > 1 cm -3 for dust scattering. Field ellipticals do not normally 
contain such large amounts of interstellar matter in a significant portion of their volume 
(Goudfrooij et al. 1994), and therefore these large values present a challenge in unifying 
elliptical hosts of quasars with their inactive counterparts. 

In Table 4 we list lower limits on the ratio of the scattered to stellar flux, calculated by 
taking into account only those positive residuals that we are confident are due to scattered 
light, i.e., confirmed by their correct orientation relative to the polarization position angle. 
This procedure only works for scattering regions that are detected on kpc scales, so that their 
orientation relative to the nucleus can be used. In several objects (e.g., SDSS J2358— 0009) 
there is a central excess of light over our best-fitting stellar profile. Although this component 
is technically resolved in all such objects (i.e., it is inconsistent with a PSF), it is impossible 
to determine if this light is due to a circumnuclear starburst or to sub-kpc scattered regions 
that are made morphologically smooth by the effects of the PSF. In principle the two can be 
distinguished by comparing their colors, especially in the UV band. However, because central 
excesses are coincident with central dust obscuration and because the quality of stellar fits 
are often poorest near the center, this method cannot be satisfactorily used for our data. 

We can also place strong lower limit on a scattered component based on the high ob- 
served polarization of the objects in our sample. The median observed polarization in the 
blue band is about 3% (Table 2). Only the scattered component is polarized, and the po- 
lazed signal is diluted by the stellar light. The polarization of the scattered light in AGNs 
has rarely been observed to be more than about 30%, and is of course limited to be <100%. 
Therefore, to produce the observed 3% polarization, at least 10% of the observed flux should 
be due to scattered light. This argument demonstrates that our photometric procedure did 
not work satisfactorily in SDSS J1106+0357 and SDSS J2358— 0009, where it failed to detect 
most of the scattered light. In SDSS J1039+6430, the high ratio of the scattered to stellar 
light (> 0.8 in the blue) is consistent with the high observed level of polarization (10% in 
the same band). 

Another method to estimate the contribution of scattered light is to include a power- 
law component in the stellar synthesis analysis of the AGN+host spectra. This method was 
employed, for example, by Kauffmann et al. (2003), who set an upper limit of a few per cent 
on the scattered light in the spectra of type II AGN hosts in their sample. This finding is 
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unsurprising given the much lower typical luminosities of the AGNs in their sample - for a 
scattering efficiency of one per cent or so, only the most luminous AGNs are expected to have 
a significant contribution of scattered light to their total flux which is otherwise dominated 
by the host galaxy and narrow emission lines. 

We used a simplified version of this type of analysis in Paper II, where we attempted 
to represent the high signal-to-noise ratio MMT spectra of several type II quasars with a 
combination of a spectrum of an old elliptical galaxy, an A-star spectrum and a power-law 
component (representing the possible scattered light). In that paper, we found that 10-100% 
of the broad-band flux was due to scattered light. The three objects that overlap between 
Paper II and this paper are SDSS J1039+6430, SDSS J1323-0159 and SDSS J1413-0142; 
the scattered-to-stellar ratio in the yellow band published in Paper II for these objects 
was estimated to be 1.0, 1.0 and 0.7, correspondingly. The spectral analysis yielded much 
higher values than those listed in Table 4, supporting our understanding that the aperture 
photometry of the positive residuals underestimates the true scattered flux (the aperture of 
the MMT spectra may have not enclosed all of the extended stellar flux of the host, but this 
effect probably cannot explain all of the difference in measured scattered-to-stellar ratios). 
Only for SDSS J 1039+6430 do the results of the two analyzes marginally agree, perhaps 
because of the lack of strong dust obscuration in the nucleus and the high quality of the 2-D 
fits. 

One of the most surprising results of our study is the great variety of shapes, sizes and 
opening angles of the scattering regions that we see in Figures 2-3. Some of them appear 
as regular, edge-brightened cones (e.g., SDSS J1323— 0159), while others appear as clumpy 
irregular patches (e.g., SDSS J1039+6430). Some scattering regions are biconical, whereas 
others are one-sided. The sample is divided almost equally between objects with large 
opening angles of scattering (30—60°) and those with jet-like, narrow regions with opening 
angles less than 10°. 

In seven galaxies, we detect kpc-scale dust extinction features. This finding is reminis- 
cent of that by Malkan et al. (1998) of the excess galactic dust in Seyfert 2 galaxies. The 
values of extinction that we find rarely exceed Ay = 1 mag, much less than what is required 
to obscure the central AGNs. (If dust is concentrated in optically thick clouds, our measured 
spatially averaged extinction is not equivalent to the probable optical depth toward the AGN 
which may in this case be much larger.) More importantly, the scattering regions seen in 
our images have a well-defined conical or bi-conical structure with opening angles that are 
much smaller than 180° - these features are hard to explain if obscuration occurs in thin 
lanes on kpc scales. Finally, in some objects (e.g., SDSS J0920+4531) the observed areas 
of extinction are not centered on the nucleus, and in two objects (SDSS J1243— 0232 and 
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SDSS J2358— 0009) we did not find any evidence for extinction in the host. Therefore, we 
conclude that this kpc-scale dust is not the primary agent of AGN obscuration, in accord 
with the standard unification picture in which the obscuration occurs on much smaller scales. 
While Malkan et al. (1998) found an interesting excess of dustiness in Seyfert 2 galaxies com- 
pared to Seyfert 1 galaxies, it is currently impossible to address the presence of extinction 
in the central parts of luminous type I quasar hosts. 

In Figure 16 we show the difference between the aperture magnitudes (Table 3) and 
the model magnitudes obtained as a result of our 2-D fitting procedure (Table 4). There 
are no apparent systematic differences between aperture and model photometry; in fact, the 
median model magnitudes and colors agree with those derived from aperture photometry 
within our claimed photometric uncertainty. It seems that, for a median galaxy, the effects 
of scattering and obscuration cancel each other out, somewhat by accident, as both effects 
become stronger toward shorter wavelengths. However, there is a large scatter in both 
luminosity differences and the color differences. If we were to use aperture photometry as a 
proxy for the luminosities and colors of the stellar component of type II quasar hosts (e.g., 
if detailed morphological information were not available), it would be accurate to no better 
than 0.3 mag. 



6. Conclusions 

In this paper, we analyze three-band ACS HST images of the host galaxies of nine 
optically-selected type II (obscured) quasars in combination with polarimetric data. Since 
optical radiation from type II quasars is heavily obscured, there is no contamination from 
the bright central source, and all host galaxies are easily detected. We further avoided the 
contribution from the narrow-line region by observing in filters that fall between strong 
emission lines. 

Type II quasars reside in elliptical hosts (6/9), in a disk+bulge host (1/9) and in dis- 
turbed disk+bulge hosts (2/9). One host can be described as undergoing a major merger, 
while up to three others show tidal debris. Type II quasar hosts are morphologically and dy- 
namically similar to the hosts of ordinary radio-quiet quasars of similar luminosities (Bahcall 
et al. 1997; Boyce et al. 1998, 1999; Hamilton et al. 2002) and of reddened quasars (Mar- 
ble et al. 2003, but see also Hutchings et al. 2003 who argue for a much higher incidence 
of interacting systems among reddened quasars). The objects presented in this paper are 
estimated to have intrinsic luminosities in the range —24 > Mb > —26 and are therefore 
of somewhat lower luminosity than the most extremely luminous unobscured radio-quiet 
quasars that seem to lie almost exclusively in massive elliptical hosts with small ellipticities 
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(Kukula et al. 1999; Dunlop et al. 2003; Floyd et al. 2004). 

The stellar light from type II quasar hosts is brighter and bluer than M* galaxies at 
the same redshift. The stellar population of the median galaxy in our sample has colors 
consistent with a 10 9 year-old starburst, while the bluest galaxy may be as young as 10 s years 
old. Several surveys have found that AGN host galaxies are significantly bluer than inactive 
elliptical galaxies. For example, based on stellar population synthesis modeling, Kauffmann 
et al. (2003) estimated that many of the type II AGN host galaxies have undergone significant 
star formation within the last (1 — 2) x 10 9 years. Our crude age estimates for the hosts of 
the high-luminosity analogs of the objects in their sample are in remarkable agreement with 
these values. Determining colors of the hosts of luminous type I AGNs is a very challenging 
task, as the already overwhelming nucleus becomes even more prominent toward shorter 
wavelengths. Studies of the hosts of medium-luminosity quasars (—23 > M B > —24) show 
that the elliptical hosts are as blue as the spiral hosts (Jahnke et al. 2004; Sanchez et al. 
2004). We find that hosts of type II quasars are still blue at intrinsic nuclear luminosities 
of —24 > Mb > —26. At yet higher nuclear luminosities, the colors and ages of quasar 
hosts are controversial. Kirhakos et al. (1999) found that even at these extreme luminosities 
the hosts of radio-loud quasars are significantly bluer than inactive galaxies of the same 
morphology and luminosity. Other authors (Kukula et al. 1999; Nolan et al. 2001) have 
argued that hosts of very luminous quasars, regardless of their radio properties, assembled 
around redshift 2.5 and have evolved passively ever since. In some of the same objects, 
Canalizo & Stockton (2000) found prominent young stellar populations and pointed out the 
importance of aperture effects in determining the ages of the stellar populations. 

Scattered light contaminates the observed images and makes it difficult to study the 
properties of the stellar populations in the host galaxies. In extreme cases, we find that 
scattered light can be the dominant source of the extended flux, especially in the UV. We 
have considered various methods to address this contamination, such as morphological anal- 
ysis (aided by polarimetric data) or stellar population synthesis analysis that includes a 
power-law template to account for a possible scattered component. For nearby well-resolved 
objects, imaging polarimetry can yield a robust determination of the extent and relative 
contribution of the scattered component (Capetti et al. 1995; Kishimoto 1999; Kishimoto et 
al. 2002a,b; see also Tadhunter et al. 2000 for imaging polarimetry of radio galaxies). At 
higher redshifts, this analysis becomes increasingly more difficult and is only possible for the 
largest scattering regions (Tran et al. 1998; see also review of polarimetric studies of radio 
galaxies by Tadhunter 2005). 

The difficulties of isolating the scattered component are expected to be more severe 
in studies of hosts of unobscured quasars where the presence of a very bright blue nuclear 
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source makes detailed morphological analysis difficult and greatly dilutes the polarization 
signal. Furthermore, since these objects are presumably seen more face-on with respect to 
the obscuration plane, there is no preferred scattering direction in the plane of the sky, 
resulting in a geometrical cancellation of the polarization from different directions. As a 
result, scattered flux can be significant without producing high polarization. In quasars with 
dust scattering (e.g., SDSS J1039+6430), the forward nature of dust scattering can make the 
contamination by the scattered light even stronger in type I quasars than in type II quasars, 
because the typical scattering angles are smaller in face-on than in edge-on objects. 

The differences between the hosts of radio-loud and radio-quiet quasars have been ad- 
dressed in many studies, but the origin and the nature of differences are still subject to 
debate. In this paper, we focused on radio-quiet type II quasars, comparing our results 
exclusively to other studies of radio-quiet quasars, unless explicitly stated otherwise. The 
radio-loud analogs of type II quasars from our sample, Narrow-Line Radio Galaxies, show 
many of the same features that we discuss in this paper, such as extended scattered light, 
significant contribution from the young stellar populations (Tadhunter et al. 2002), and a 
range of morphologies and dynamical states (de Kofi et al. 1996). A careful statistical com- 
parison of large samples of radio-quiet and radio-loud type II AGNs is needed to address the 
spectacular difference in their radio output. 

In the Introduction, we postulated that the hosts of type II AGNs form a representative 
sample of all AGN hosts. Statistical differences between type II and type I hosts have been 
claimed in the literature in the past, suggesting that some connection is present between 
the parsec-scale nuclear region and the kpc-scale host (Antonucci 1989). Such comparison 
becomes increasingly more complicated at higher luminosities, but so far we have not seen 
any dramatic differences between the morphologies, dynamical states, luminosities and ages 
of the stellar populations of the hosts of type II vs type I AGNs. 

From spectropolarimetric observations presented in Paper II, we know that at least 
some, if not all, of our type II quasars are seen as broad-line AGNs along some directions. 
Therefore, our findings apply to at least some type I quasars. The question remains open 
as to whether all type I quasars are obscured along some lines of sight or whether there 
exists a population of truly 'naked' type I quasars that are seen without circumnuclear 
dust obscuration along any direction. The calculation of the type I/type II ratio at high 
luminosities meets serious statistical difficulties, because there is at present no robust method 
of AGN selection that is not biased toward a specific type. 
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A. Broad-band polarimetry 

Broad-band polarimetric measurements for 11 objects from Paper I are presented in 
Table 5. Observations were conducted using a V filter (4800 - 6000 A). The measured values 
are Q, U and a(P). The derived values are tan 29 = U/Q (where 9 is in degrees East 
of North), P = ^Q 2 + U 2 ) a{0) = 28.65° x a(P)/P, and q = V ' P 2 - a(P) 2 (only for 
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the objects with P > cr(P)). The polarization position angle is essentially undefined if 
P < 1.5 x <j(P) (Wardle k Kronberg 1974). The value P is positive definite and therefore in 
the presence of measurement errors is an overestimate of the true polarization level q. The 
latter is estimated using the approximation by Wardle k Kronberg (1974). The observational 
setup is described in Section 2.3. With the exception of SDSS J0056+0032, polarization is 
detected in all objects at less than the 2a level. 
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Fig. 1. — SDSS spectra of the nine objects that are the subject of this paper. The resolution 
of the original spectra is A/ A A 2000; the spectra have been smoothed by 5 pixels for 
display purposes. Overlaid on the spectra are transmittance curves of the regular ACS 
filters (dotted lines) and ramp ACS filters (dashed lines, schematic representation) used in 
our imaging program, as summarized in Table 1. The filters were chosen to avoid strong 
emission lines. Transmittance curves and widths of ramp filters were taken from the ACS 
instrument website, http://acs.pha.jhu.edu/instrument/. 
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Fig. 2. — Color-composite images of type II quasars produced using the algorithm by Lupton 
et al. (2004). In all images except SDSS 1323—0159 (top left) the same parameters of the 
color-combining code were used (group 1 objects, see Table 1), so the colors of the objects 
can be directly visually compared. Object identification and the total horizontal size of 
the images are given above each frame. A cartoon of the presumed scattering regions and 
polarization position angle (in image coordinates), as well as the orientation on the sky 
are shown to the right of each image. The polarization position angle 9 is marked with 
a solid black line, and grey lines mark the lcr confidence limits if the uncertainty in the 
angle is significant. The brightest UV regions tentatively identified as scattered light are 
schematically shown as grey wedges. 
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Fig. 3. — Color-composite images of type II quasars produced using the algorithm by Lupton 
et al. (2004). In all three images, the same parameters of the color-combining code were 
used (group 2 objects, see Table 1), although the positions of the ramp filters used for 
SDSS J1243— 0232 were slightly different from those used for the other two objects (Figure 
1 and Table 1). The UV filter used for the three objects in this image (FR459M) is about 
an order of magnitude less sensitive than F435W used in Figure 2, so the scattering regions 
are harder to identify, and shot noise can be seen in the form of blue spots on top of a dark 
background. Images are labeled as in Figure 2. We do not have polarimetric measurements 
of SDSS J1243-0232 and SDSS J1301-0058. 
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Fig. 4. — Results of spectropolarimetry conducted using SPOL on MMT for 
SDSS J0920+4531 (left) and SDSS Jl 106+0357 (right). From top to bottom: (a) total 
optical flux density (focusing on the continuum); (b) polarized flux density; (c) polarization 
fraction; (d) polarization position angle in degrees east of north. Flux densities in panels (a) 
and (b) are in units of 10~ 17 erg sec -1 cm -2 A -1 . All measurements are corrected for Galactic 
extinction and atmospheric absorption bands. Polarization measurements have been heavily 
binned in wavelength. The horizontal error bars reflect the bin size, whereas the vertical 
error bars are la uncertainties computed from the scatter of values within each bin. 
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Fig. 5. — 1-D brightness profiles of the extended emission, regardless of its origin, for all 
objects in all bands. For each object, log J(r) is plotted in units of apparent AB mag/arcsec 2 
as a function of the semi-major axes of the isophotes r and as a function of r 1//4 . Yellow 
band brightness profiles are shown in yellow, blue band in blue and UV in purple. The error 
bars reflect the rms scatter of the values of the surface brightness on the fitted isophote. 
The routine fails in the outer parts of SDSS J1301— 0058 because of the presence of a nearby 
companion, so the profiles are shown only for the inner 4 kpc. The UV filters used for 
SDSS J1106+0357, SDSS J1243-0232 and SDSS J1301-0058 were significantly less sensitive 
than those used for other objects, so the UV 1-D profiles of these objects cover only the well- 
detected inner parts. 
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Fig. 6. — Examples of fitting elliptical isophotes for three objects. From left to right: (a) the 
image, (b) the model reconstructed from the best-fit elliptical isophotes, (c) residuals. As 
seen in the top panel, residuals in the yellow are 1% or less except for a few central pixels, 
and are typically much lower than those in the UV (bottom two objects) where they reach 
tens of per cent. All black-and-white images in this paper are negatives: darker areas mean 
more flux per pixel. Orientation of these and subsequent images is the same as in Figures 
2-3. 
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Fig. 7. — Subtracting a 2-D de Vaucouleurs profile from the image of SDSS JO 123+0044. 
From top to bottom: yellow, blue and UV bands, on the same scale (the size of each frame 
is 70 pixels, or 18.4 kpc). From left to right: the image, the best-fitting model, and the 
residuals. Residuals become stronger toward shorter wavelengths, and the uncertainty of 
the model increases. In particular, the differences in the position angles and ellipticities 
between the bands are within the model uncertainties. Both positive and negative residuals 
are present; the latter appear as solid white. While negative residuals are probably due to 
dust obscuration, the nature of the positive residuals is less certain. 
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Fig. 8. — Faint extended emission associated with SDSS J0123+0044 in the UV band (the 
tail on the left-hand side of the image; the feature is not detected in the yellow and blue 
bands). The image has been smoothed with an 8x8 boxcar filter for display purposes. The 
horizontal scale of the frame is 35.8 kpc. 



Fig. 9. — Subtracting a de Vaucouleurs profile from the image of SDSS J1039+6430. The 
frames, which are 82 pixels (21.6 kpc) wide, are in the same order as in Figure 7. Almost all 
residuals are positive and are due to scattered light. 
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Fig. 10. — Subtracting a sum of a de Vaucouleurs profile and exponential profile from the 
image of SDSS J1106+0357 in the yellow band. The frames, which are 254 pixels (48.5 
kpc) wide, show the original image and the residuals after subtraction. Both frames have 
been smoothed with an 8x8 boxcar to highlight the low surface brightness features. A 
faint extended arm connecting the main galaxy to an extended companion to the north-east 
(image left) of the nucleus can be seen, especially in the residual image. 



-41 - 




Fig. 11. — Subtracting a de Vaucouleurs profile from the image of SDSS J1243— 0232. The 
frames, which are 130 pixels (27.1 kpc) wide, are in the same order as in Figure 7. There 
is excess emission in the center relative to the best-fit de Vaucouleurs profile; this excess is 
also apparent in the 1-D profile of this object (Figure 5). This object is one of the sources 
which was imaged using the low-sensitivity ramp filter FR459M (Table 1, group 2), and the 
extended UV emission is just barely detected. 
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Fig. 12. — Subtracting a de Vaucouleurs profile from the image of SDSS J1323— 0159. The 
frames, which are 55 pixels (13.3 kpc) wide, are in the same order as in Figure 7. The 
scattering cones (X-shaped structures) are apparent in the residuals. The central flux is 
strongly suppressed relative to the best-fitting de Vaucouleurs profile, in agreement with the 
results of the 1-D analysis (Figure 5). This is indicative of a kpc scale obscuration in this 
object, which is also suggested by the color-composite image (Figure 2). 
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Fig. 13. — Faint extended emission associated with SDSS J2358— 0009 in the yellow band. 
The image has been smoothed with an 8x8 boxcar. The horizontal scale of the frame is 82.4 
kpc. The two black blobs are the two galaxies visible in Figure 2. 
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Fig. 14. — (a) Distribution of absolute luminosities of type II quasar hosts relative to the 
M* values in the yellow band (solid; median value — 0.3loi mag; bootstrap error on the 
median was used) and in the blue band (dotted; median value —0.7^0% mag). Values < 
mean brighter than M*. (b) Distribution of colors of type II quasar hosts relative to the M* 
values: blue— yellow colors (solid; median value —0.4 mag) and UV— yellow colors (dotted; 
median value —0.2 mag). Relative colors < mean bluer colors than those of M* galaxies. 
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Fig. 15. — Nuclear vs host luminosities of type II quasars (open circles with error bars) and 
of unobscured quasars (blue: radio-quiet, red: radio-loud). Different symbols correspond to 
objects from different studies (filled squares: Hamilton et al. 2002, open triangles: Bahcall 
et al. 1997, crosses: Floyd et al. 2004, open squares: Jahnke et al. 2004). All objects from 
the sample by Bahcall et al. (1997) were re-analyzed by Hamilton et al. (2002); for a given 
object, the results of the two analyzes are connected with a solid line. All magnitudes from 
the literature were corrected to the h = 0.7, Q m = 0.3, = 0.7 cosmology. Dotted lines 
are the contours of constant total (nucleus +host) flux. Intrinsic nuclear luminosities of type 
II quasars are estimated using the [OIII]5007— My correlation that exists for unobscured 
quasars (Paper I). 
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Fig. 16. — (a) Distribution of model magnitudes of type II quasar hosts relative to the 
aperture magnitudes in the yellow band (solid; median value —0.1 mag), in the blue band 
(dotted; median value —0.1 mag) and in the UV band (grey histogram; median value 0.1 
mag). Values < mean model magnitudes brighter than aperture magnitudes, (b) Distri- 
bution of model colors of type II quasar hosts relative to the aperture values: blue— yellow 
colors (solid; median value 0.0 mag) and UV— yellow colors (dotted; median value 0.1 mag). 
Relative colors > mean model colors are redder than aperture colors. 



Table 1. Summary of HST observations 



Object name 


redshift, 


log 


filters 


rest- frame eff. A, A 






(J2000) 


z 


L[OIII]/L 


yellow, blue, UV 


yellow 


blue 


UV 


polarimetry 


group 


SDSS J012341.47+004435.9 


0.399 


9.13 


FR914M, F550M, F435W 


5700 


3989 


3086 


S02.3 


1 


SDSS J092014.11+453157.3 


0.402 


9.04 


FR914M, F550M, F435W 


5688 


3981 


3079 


MMT 


1 


SDSS J103951.49+643004.2 


0.402 


9.41 


FR914M, F550M, F435W 


5688 


3981 


3079 


MMT (Paper II) 


1 


SDSS J110621.96+035747.1 


0.242 


9.13 


FR647M, F550M, FR459M 


5804 


4494 


3124 


MMT 


2 


SDSS J124337.34-023200.2 


0.281 


9.02 


FR647M, F550M, FR459M 


5789 


4357 


3118 




2 


SDSS J130128.76-005804.3 


0.246 


9.25 


FR647M, F550M, FR459M 


5786 


4479 


3114 




2 


SDSS J132323.33-015941.9 


0.350 


9.19 


F775W, F550M, F435W 


5699 


4134 


3198 


MMT (Paper II) 




SDSS J141315.31-014221.0 


0.380 


9.25 


FR914M, F550M, F435W 


5779 


4044 


3129 


MMT (Paper II) 


1 


SDSS J235818.87-000919.5 


0.402 


9.32 


FR914M, F550M, F435W 


5688 


3981 


3079 


S02.3 


1 



Note. - - J2000 coordinates, redshifts and [OIII]5007 luminosities are from the SDSS data (Paper I). For five objects 
spectropolarimetry data from MMT are available; for three of these the data were published in Paper II, and for the remaining 
two the data are presented in this paper. For two additional objects broad-band polarimetry obtained using the Steward 
2.3m telescope (S02.3) is available. In the last column, we grouped the objects for display purposes based on similarity of 
the filters used in the observations. Objects from the same group have the same parameters of the color-composite (Lupton 
et al. 2004), and therefore can be visually compared with one another (Figures 2-3). 



Table 2. Summary of polarization properties 



Object 


UV 


blue 


yellow 


mean 9 


name 


polarization polarization polarization de 


g. E of N 


SDSS J0123+0044 




2.7±2.2 




8± 25 


SDSS J0920+4531 


4.7±0.9 


2.5±0.4 


0.7±0.6 


171±2 


SDSS J1039+6430 


16.6±0.3 


10.2±0.2 


8.0±0.3 


109±1 


SDSS J1106+0357 


3.1±0.5 


1.0±0.3 


0.1±0.2 


146±3 


SDSS J1323-0159 


5.0±1.0 


3.7±0.5 


2.8±0.5 


104±3 


SDSS J1413-0142 


4.1±1.0 


2.9±0.4 


1.2±0.6 


146±8 


SDSS J2358-0009 




2.8±1.8 




154±15 



Note. - 'UV, 'blue' and 'yellow' polarizations (given in %) are 
calculated as close to the ACS filter coverage as possible. We list the 
mean polarization position angle 9 and its la error (calculated over the 
entire available polarimetric wavelength coverage). 



Table 3. Aperture photometry of type II quasars 





Object 


absolute AB mag 


nitudes 




name 


vp 1 low 

yciiu vv 


"hlnp 


UV 


oDbo 


JU123+UU44 


-21.8 


-20.7 


-19.9 


SDSS 


J0920+4531 


-22.4 


-21.6 


-20.7 


SDSS 


J1039+6430 


-22.0 


-21.6 


-21.0 


SDSS 


J1106+0357 


-22.7 


-22.0 


-19.6 


SDSS 


J1243-0232 


-22.6 


-21.9 


-19.9 


SDSS 


J1301-0058 


-21.7 


-21.1 


-19.3 


SDSS 


J1323-0159 


-21.4 


-20.7 


-19.8 


SDSS 


J1413-0142 


-21.6 


-20.9 


-19.9 


SDSS 


J2358-0009 


-22.5 


-21.8 


-20.9 



Note. - Magnitudes are given at effective 
wavelengths listed in Table 1. 



Table 4. Structural parameters and luminosities of type II quasar host galaxies 



Object 


Sersic 




Rd 


or R e , kpc 




mo 


del AB ma; 


mitudes of the stellar component 


scattered-to- 


name 


index 


ell. 


yellow 


blue 


UV 


yellow 


(yellow*) 


blue 


(blue*) 


UV 


(UV*) 


stellar 


SDSS J0123+0044 


4 


0.45 


3.7 


6.0 


4.9 


-21.9 


(-21.7) 


-20.8 


(-20.6) 


-19.7 


(-19.9) 


> 0.1 


SDSS J0920+4531 


4+1 


0.6 








-22.4 


(-21.7) 


-21.7 


(-20.6) 


-20.7 






SDSS J1039+6430 


4 


0.1 


1.4 


2.3 


1.9 


-21.4 


(-21 7) 


-21.1 


(-20.6) 


-20.1 


(-19.9) 


> 0.8 


SDSS J1106+0357 


1 


0.45 


11.0 


14.3 




-22.9 


(-21.5) 


-22.5 


(-20.6) 


<-19.2 


(-19.2) 


>0.01 




4 


0.4 


5.5 


3.8 


1.6 


total 




total 




total 






SDSS J1243-0232 


4 


0.1 


6.5 


7.0 


8.1 


-22.5 


(-21.6) 


-21.7 


(-20.6) 


-19.9 


(-19.4) 




SDSS J1301-0058 


4 


0.2 


10.3 


14.4 




-21.8 


(-21.5) 


-21.3 


(-20.6) 




(-19.2) 




companion 


4+1 


0.5 












-19.2 










SDSS J1323-0159 


4 


0.35 


1.9-2.6 


1.2-3.2 


0.6-1.8 


-21.5 


(-21.7) 


-20.9 


(-20.6) 


-20.1 


(-19.7) 


>0.08 


SDSS J1413-0142 


4+1 


0.45 








-21.5 


(-21.7) 


-20.5 


(-20.6) 


-19.8 


(-19.8) 


>0.09 


SDSS J2358-0009 


4 





8.6 


11.1 




-22.7 


(-21.7) 


-21.9 


(-20.6) 


-20.8 


(-19.9) 


>0.005 Or 


companion 


4+1 











-21.6 




-20.7 




-19.0 




o 



Note. — Parameters of the best 2-D model fits (de Vaucouleurs if the Sersic index is 4, exponential if it is 1). The column 'ell.' lists our best-fit 
ellipticities in the blue band. Typical uncertainties in the characteristic radius are 20% or more. Absolute model magnitudes are given at the 
effective rest- frame wavelengths listed in Table 1. The values with the * symbol are absolute magnitudes of the M* galaxies at the same redshift 
and the same effective wavelength (Section 5.1). Model magnitude uncertainties are 0.1 mag except for SDSS J0920+4531, SDSS J1323— 0159 (UV) 
and SDSS J1413— 0142. In these cases, we estimate the uncertainties to be about 0.2-0.3 mag. In SDSS 1301—0058, only the central excess emission 
of unknown nature was detected in the UV band, and in SDSS Jl 106+0357 only the bulge component was detected in the UV. In the last column, 
we list lower limits on the ratio of the scattered to stellar fluxes in the blue band, based on the aperture photometry of the positive residuals 
confidently identified as scattered regions. In SDSS J0920+4531, bright positive residuals were detected, but their orientation does not agree with 
the polarization position angle. Polarization data were not available for SDSS J1301— 0058 and SDSS J1243— 0232, so none of the residuals could 
be identified as scattered light. 
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Table 5. Appendix: polarization measurements of 11 type II AGNs from Paper I 





J2000 








log 


Q 


U 


P 




Q 


a(P) 


e 


a(6) 




coordinates 






z 


L[OIII]/L 


% 


% 




Vo 




% 


% 


deg. 


deg. 


SDSS 


J002531.46-104022. 


.2 





,303 


8.73 


-0.95 


-0.80 


1 


.24 


1 


,07 


0.62 


110 


14 


SDSS 


J002827.78-004218 


.8 


0. 


.418 


8.75 


-2.17 


-1.67 


2 


.74 


1 


.47 


2.31 






SDSS 


J005009.81-003900. 


.6 





,729 


9.94 


2.93 


1.10 


3 


.13 


2 


.45 


1.94 


10 


18 


SDSS 


J005621. 72+003235. 


,8 


0. 


,484 


9.45 


0.72 


-8.66 


8 


.69 


8 


.31 


2.56 


137 


8 


SDSS 


J012341. 47+004435. 


.9 


0. 


.399 


9.13 


3.33 


0.97 


3 


.47 


2 


.67 


2.21 


8 


18 


SDSS 


J014237.49+144117. 


.9 


0. 


,389 


8.76 


-0.64 


-0.16 





.66 






0.76 






SDSS 


J015911.66+143922. 


,5 


0. 


.319 


8.56 


-0.02 


-2.62 


2 


.62 






3.64 






SDSS 


J021758. 19-001302 


.7 





.344 


8.75 


0.53 


-1.43 


1 


.53 


1 


.18 


0.97 


145 


18 


SDSS 


J225102.40-000459 


.9 


0. 


.550 


9.13 


-2.87 


-1.65 


3 


.31 


2. 


.69 


1.93 


105 


17 


SDSS 


J235818.87-000919. 


.5 


0. 


.402 


9.32 


2.01 


-2.61 


3 


.30 


2. 


.79 


1.75 


154 


15 


SDSS 


J235831. 16-002226. 


.5 


0. 


,628 


9.96 


4.75 


1.19 


4 


.90 


4 


.23 


2.47 


7 


14 



Note. — Broad-band polarization measurements for 11 objects from Paper I. Coordinates, 
redshifts and [OIII] luminosities (quoted as logL[OIII]/L Q ) are from Paper I. Polarization 
values are quoted as P (uncorrected for statistical bias) and q (corrected for statistical bias). 
is the polarization position angle, measured East of North. 



